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Thermal and Fragmentation Properties of Star-forming Clonds 
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ABSTRACT 

The thermal and chemical evolution of star-forming clouds is studied for different 
gas metallicities, Z, using the model of Omukai (2000), updated to include deuterium 
chemistry and the effects of cosmic microwave background (CMB) radiation. HD-line 
cooling dominates the thermal balance of clouds when Z ~ 10“^ — 10“^Zq and density 
10^cm“^. Early on, CMB radiation prevents the gas temperature to fall below Tcmb, 
although this hardly alters the cloud thermal evolution in low-metallicity gas. From 
the derived temperature evolution, we assess cloud/core fragmentation as a function 
of metallicity from linear perturbation theory, which requires that the core elongation 
S = {b — a)/a > <Snl ~ 1, where a (b) is the short (long) core axis length. The fragment 
mass is given by the thermal Jeans mass ai £ = Znl- Given these assumptions and 
the initial (gaussian) distribution of £ we compute the fragment mass distribution as 
a function of metallicity. We find that: (i) For Z = 0, all fragments are very massive, 
< 10^ Mq, consistently with previous studies; (ii) for Z > 10“®Zq a few clumps go 
through an additional high density (> 10^°cm“^) fragmentation phase driven by dust¬ 
cooling, leading to low-mass fragments; (iii) The mass fraction in low-mass fragments 
is initially very small, but at Z ~ 10' -5Zo it becomes dominant and continues to grow 
as Z is increased; (iv) as a result of the two fragmentation modes, a bimodal mass 
distribution emerges in 0.01 < Z/Zq < 0.1. (v) For > O.IZq, the two peaks merge into 
a singly-peaked mass function which might be regarded as the precursor of the ordinary 
Salpeter-like IMF. 

Subject headings: cosmology: theory — galaxies: formation — stars: formation 


1. Introduction 

The presence of a massive stellar population in the early universe has been frequently invoked 
to account for some recent observations. Most notably the large electron scattering optical depth 
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of Te = 0.17+0;[;® detected by the WMAP satellite, corresponding to the reionization of the inter- 
galactic medium at redshifts < 10 (Spergel et ah 2003; Kogut et al. 2003), is interpreted with 
photoionization by massive population III (Pop III) stars. In addition, the amplitude and spectral 
break found at 1 /xm in the infrared background radiation is explained by radiation from Pop III 
stars at z > 9 (Santos, Bromm, & Kamionkowski 2002; Salvaterra & Ferrara 2003). 

Even before those evidences, theorists claimed that the first stars in the universe should be 
very massive, typically 100 — IOOOMq (e.g., Abel, Bryan, & Norman 2000, 2002; Omukai & Palla 
2001, 2003). Among them, Bromm et al. (1999, 2002) investigated the collapse and fragmentation 
of primordial-gas clouds with 3-dimensional numerical hydrodynamics. They discovered that frag¬ 
mentation occurs at density ~ 10^cm“^. As a result of this, dense cores with mass-scale ~ IOOOMq 
are produced, which do not fragment into smaller sub-clumps after that. This fragmentation history 
can be interpreted from the view point of thermal evolution of primordial gas. In low-temperature 
(< lO^K) primordial gas, the only important cooling agent is rotational/vibrational line emission 
by molecular hydrogen. Once a sufficient amount of H 2 for dynamical collapse is formed, resultant 
cooling causes a temperature decrease with increasing density. The collapse in this phase is so 
rapid that density perturbations have not enough time to grow and the clouds do not fragment. 
However, at ~ 10^cm“^ the rotational levels of H 2 reach local thermodynamic equilibrium, and the 
cooling becomes inefficient at higher densities. Now the temperature begins to increase with in¬ 
creasing density. The transition from temperature decrease to increase causes the clouds to remain 
for long time intervals (or to “loiter”) at ~ 10^cm“^ at 200K, the minimum temperature needed 
for excitation of H 2 levels (the hrst excited level locates at ~ 600K). At this moment, as there is 
sufficient time for fluctuations to grow, fragmentation occurs. The fragmentation mass-scale is set 
by the thermal Jeans mass at this time, i.e., ~ IOOOMq, using the number density ~ 10^cm“^ and 
temperature ~ 200K. Those massive dense cores are expected to grow eventually to massive stars 
because even though protostars are initially small, rapid accretion onto the formed protostars is 
not halted until at least it grows to ~ IOOMq without dust grains (Omukai & Palla 2003; Tan & 
McKee 2004; Bromm &: Loeb 2004). 

Along the same lines, the characteristic mass-scale (0.1 — IMq) of present-day Galactic star 
formation may be derived (Larson 2004). In the high-mass regime (< Mq), the number of field 
stars decreases with increasing stellar mass following the Salpeter IMF, while the slope of the IMF 
flattens in the sub-solar mass regime and eventually decreases for stars below about O.IMq. The 
precise origin of this mass-scale is still obscure, but Larson (2004) interpreted it as corresponding 
to the Jeans mass at the end of rapid cooling phase caused by dust thermal emission. 

From the arguments above, a close relation between the thermal evolution and fragmentation 
is evident. Fragmentation occurs at the end of temperature-decreasing phase, or in other words, 
when the isothermality breaks, namely the effective ratio of specific heat 7 = dlogp/dlogp > 1 ; the 
typical fragmentation mass-scale is given by the thermal Jeans mass at that epoch. 

Significance of the role played by the massive stellar population in early universe crucially 
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depends on the duration of massive-star formation epoch and the amount of massive stars formed. 
As metal enrichment proceeds in interstellar medium (ISM), it might affect the stellar mass-scale in 
two different ways. First, pre-stellar thermal evolution is altered, and consequently the typical mass- 
scale of fragmentation changes with metallicity. Second, in sufficiently high dust-to-gas ratio cases, 
i.e. high-metallicity environment, the accretion onto massive protostars is inhibited by radiation 
force onto dust grains resulting in the scarcity of massive stars. As for the former effect, even 
with metallicities as small as 10 ’ the thermal evolution is altered significantly, and low-mass 

fragments are likely to be produced (Schneider et al. 2002 ; 2003). As for the latter, a metallicity as 
high as ~ 10~^Zq is necessary to affect the accretion process by radiation force onto dust grains, if 
the dust opacity for a given metallicity remains the same as in local ISM (Omukai & Palla 2003). 
Since the former effect becomes important earlier in the history of metal-enrichment, the transition 
from massive to low-mass star formation mode is probably caused by the fragmentation process 
rather than by halting the accretion. 

Based on the thermal evolution obtained by Omukai (2000; Paper I), Schneider et al. (2002, 
2003) discussed fragmentation properties for low-metallicity clouds and pointed out that even within 
the metallicity range 10 “^^^.^0 low-mass fragmentation is possible. However, the thermal evolution 
of those low-metallicity clouds is rather complicated. The molecular cooling phase is followed by 
temperature increase owing to H 2 formation, and finally by dust-induced temperature drop before 
becoming optically thick. That is, there are a few up-and-downs in temperature evolution and the 
epoch of isothermality break is not unique. Therefore, it is difficult to determine whether fragmen¬ 
tation indeed occurs. In particular, after the first fragmentation in molecular-cooling phase, the 
cloud starts run-away collapse, during which density fluctuations might be erased by pressure force. 
If so, with little elongation of the core, the fragmentation at higher density would be prohibited. 
Therefore, whether fragmentation occurs in the dust-cooling phase is still unclear. 

To clarify this point, here we follow the evolution of elongation of dense cores during the 
run-away collapse and discuss when the final fragmentation occurs. Our results confirm that frag¬ 
mentation does occur in the dust-cooling phase and sub-solar mass fragments are indeed produced. 
In some cases, dust-induced fragmentation occurs even in metallicity as low as 10~^Zq. We have 
improved the previous model (Paper I) by including deuterium chemistry. Although HD cooling is 
important neither in the metal-free case nor in metal-rich cases, HD is found to be the dominant 
coolant in cases with slight (~ 1 O“^.^ 0 ) metallicity for some density ranges. Finally, we have con¬ 
structed reduced chemical models that can reproduce the thermal evolution in an approximate way 
and might be used to include the chemical network of metals in future challenging simulations. 

The paper is organized as follows. In §2, our model for star-forming clouds is briefly described 
with emphasis on modification from the previous one. In §3, the calculated thermal and chemical 
evolution is presented. Finally in §4, we summarize and discuss the results with some speculations 
on the evolution of IMF. In Appendix, by applying a simple model for fragmentation to the thermal 
evolution derived in §3, the fragmentation mass-scales are discussed. 
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2. Model for Pre-stellar Clouds 


Our analysis in this paper is based on the model developed in paper I. Prestellar clouds 
are treated by a one-zone model where involved physical variables are regarded as those at the 
center. Thermal and chemical variables at the center, e.g., temperature T = T(nH), concentration 
of chemical species i y{i) = y(i)(nH)^ etc., are followed as a function of the number density of 
hydrogen nuclei nn- While detailed radiative and chemical processes are included in this model, 
the dynamics of the clouds is largely simplified. 


The temperature evolution is followed by solving the energy equation; 


de d 1 

It ^ ~ 


where the pressure 


P = 


pkT 

pmn 


and the specific thermal energy is written as 


( 1 ) 

( 2 ) 


1 kT 
7ad - 1 pmn 


( 3 ) 


by using the adiabatic exponent 7 ad, the temperature T, and mean molecular weight p,. In equation 
(1), the first term on the right hand side is the compressional heating, and the second term Anet is 
the net cooling due to the radiative processes and chemical reactions. The net cooling rate consists 
of three components: 

■^net — -^line T -^cont T -^chem) (d) 


i.e., the cooling rate owing to line radiation Anne, to continuum radiation Acont, and to chemical 
reactions Achem- The last term can be negative for chemical heating rather than cooling. In our 
treatment, as in equation (3), the specific thermal energy does not include chemical binding energy. 
Variations in chemical binding energy are regarded as cooling and heating processes and included 
in the net cooling rate Anet- As the line radiation process, we include molecular line emission of H 2 , 
HD, OH and H 2 O, and atomic fine-structure line emission of [Cl], [CH], [01]. The H 2 collisional 
transition rates of low-lying rotational levels (J < 5) are updated by following Galli & Palla (1998) 
from paper I where Hollenbach & McKee (1979)’s coefficients were used. This modification lowers 
the H 2 cooling rate below 300K and alters the behavior around nn ~ 10^cm“^ for low-metallicity 
clouds. The HD parameters are from Galli & Palla (1998)’s compilation. Other parameters are 
the same as in paper I. Continuum cooling consists of that by dust thermal emission and by H 2 
collision-induced emission. In optically thick cases, cooling rate is reduced by multiplying by the 
escape probability of photons. Unless otherwise stated, the treatment of these processes are the 
same as in paper I. 


^The concentration y{i) of the chemical species i is defined by y{i) = n{i)/n-H, where n(i) is the number density 


of the species i. 
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External radiation is not included in paper I. Cosmic microwave background radiation (CMB) 
in early universe, whose temperature is as high as several tens K, can significantly influence the 
thermal evolution of star-forming clouds. Although there should exist additional radiation field 
from stars or quasars, the amount and spectrum are quite uncertain. Therefore, in this paper, we 
consider only CMB as an external radiation source. CMB enters in our model in two ways; first 
through the level population of atomic/molecular lines, and second through the heating of dust 
grains. These effects prevent the gas and dust temperature to fall below the radiation temperature. 
The former effect is included in calculating level population of atomic and molecular species as in 
Omukai(2001). The latter effect is included in determining the grain temperature by the energy 
balance equation of grains: 

du'T'gj.Kgr = Agag^dust T (5) 

where Agas-dust is energy transfer rate from gas to dust per unit mass, which is given by equation 
(19) of paper I The left hand side of equation (5) is the cooling rate of grains and the right hand 
side is heating rate by collisions with gas particles (the first term) and by CMB heating (the second 
term). Here we use the approximation that the Planck mean opacity for absorption is the same 
as that for emission. This is a good approximation, being exact when Tgr = Trad) i-e. when CMB 
effect is significant. 

In paper I, we have solved H, He, C, and O chemical network among 50 species and 478 
reactions. In this work, we consider the additional contributions of deuterium chemistry and 
HD cooling. The deuterium chemical network and the reaction rates are taken from Table 1 of 
Nakamura & Umemura (2002), which include 18 reactions among 5 species (D, D"*", D“, HD, 
HD"*"). Deuterium abundance is i/d = 3 x 10“^. For gas metallicity equal to that of the solar 
neighborhood .^locab gas-phase elemental abundances of carbon and oxygen are yc,gas = 0.927 x 
10“*^ and yo,gas = 3.568 x 10“^. The mass fraction of grains is 0.939 x 10“^ below the water-ice 
evaporation temperature, and it decreases for higher temperatures as, following the grain model by 
Pollack et al. (1994), each grain component evaporates. For lower gas metallicities, these values 
are reduced proportionally. We denote the relative abundances relative to the solar neighborhood 
values by [Z/H]= logio(^/^iocai)- 

Next, we describe our prescription for the dynamics. In our model, the only two quantities 
related to the dynamics are the collapse timescale and the optical depth. Since we are interested in 
the evolution during the run-away collapse phase, we assume that the dynamics is described by the 
Larson-Penston-type self-similar solution, which reproduces the dynamics of the run-away collapse 
phase very well (Larson 1969; Pension 1969).^ Given the collapse timescale tcoi of the self-similar 

^This equation contains a typographical error. The correct numerical factor is 1.1 x 10~® instead of 1.1 x 10® 

®The solution was found by Larson (1969) and Pension (1969) for the isothermal case and later extended to the 
polytropic cases by Yahil(1983). 
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solution, the evolution of the central density is followed by 


dp 

dt 


P 

tco\ 


The collapse timescale of the self-similar solutions is given by 

1 

^col — j ^col,0 

V-L - / 

where / is the ratio between the pressure gradient force and gravity at the center and 


( 6 ) 

(7) 


I Svr 

= V 3^ 

is the collapse timescale for pressure-free collapse. We approximate the ratio of pressure to gravity 
/ in the central region of the self-similar solutions by the following formula; 


/ 


0 7 < 0.83 

< 0.6 -h 2 . 5(7 - 1) - 6 . 0(7 - 1)^ 0.83 < 7 < 1 

1 . 0 -h 0 . 2 ( 7 -4/3)- 2 . 9 ( 7 -4/3)2 7>1- 


(9) 


Error in / is less than 0.05. At 7 = 4/3, / reaches unity and tcoi diverges. This indicates a 
break-down of our assumption that hydrodynamical evolution is approximated by the self-similar 
solution. In reality, a hydrostatic core is formed at this point. Thereafter the core contracts as its 
mass increases owing to accretion of the envelope. To mimic this evolution, we set an arbitrary 
upper limit of / = 0.95. The value of this upper limit does not affects the thermal evolution 
significantly unless when it is very close to unity, say / > 0.99. Early in the evolution, before 
the clouds experience the first fragmentation, the dynamics has not yet converged to the run-away 
collapse. In this regime, we assumed free-fall collapse, i.e., / = 0. 


The density distribution of the self-similar solution, which is needed to estimate the optical 
depth, is as follows; inside the central region, whose length scale is roughly given by thermal Jeans 
length 

Aj = c^^TxjGp, (10) 

where Cg is the sound speed, the density is approximately spatially constant, while in the envelope 
(^ ^ -^j)) it decreases outward as p oc r“ 2 /( 2 - 7 ). Since the column density in the envelope is at 
most the same as that of the central region and, in addition, the large Doppler shift owing to bulk 
velocity reduces the optical depth of the lines, we consider only the shielding of photons by material 
in the central region. Namely, if the cloud becomes optically thick, the optical depth is evaluated 
by the relation 

Ty = (II) 

where Ky is the opacity at the center. 

We neglect dark matter gravity in computing the collapse timescale. Although it might be 
important for the low-density regime, at sufficiently high densities, where fragmentation is expected 
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to occur, the baryons have already contracted significantly owing to energy dissipation and dominate 
the gravity. Indeed, fragmentation is caused by the self-gravity of the baryons and the contribution 
of the dark matter does not alter the dynamics. 

We initialize the model with an arbitrary set of initial conditions for the density and temper¬ 
ature of the gas. In particular, we assume an initial number density of nn = 0.1cm“^ and, at this 
time, the temperature is set to be T = 300K from Bromm et al. (2002)’s result (see their Figure 10 
c). Indeed, the temperature evolution converges immediately to a certain trajectory, and variations 
in the initial conditions affect only the short transient phase of the initial evolution (e.g., Omukai 
2000). For this reason, in this paper we do not consider different initial conditions and we refer the 
interested reader to, e.g., Palla et al. (1983) and Omukai (2000) for a thorough discussion of this 
issue. The initial ionization degree and H2 concentration are assumed to be the values for uniform 
matter in the post-recombination universe: y{e) = 10“^, y{H 2 ) = 10“®. Initially, carbon is fully 
ionized and Helium and Oxygen are neutral. 


3. Results 

In this section, we first present the thermal and chemical evolution for different initial values 
of the gas metallicity (§3.1); then, a reduced chemical model for low-metallicity clouds is presented 
(§3.2). Using the derived effective equation of state, we discuss the fragmentation properties of 
prestellar clouds and estimate the typical mass-scales of fragmentation in Appendix. 


3.1. Thermal and Chemical Evolution 

3.1.1. Results for the Fiducial Cases 

The temperature evolution for prestellar clouds for different metallicities is shown in Figure 
1. The cases with metallicities [Z/H\ = —00 (i.e., Z=0 metal-free case), -5, -3 and -1 {[Z/H] = 
—6, —4, —2 and 0) are indicated by solid (dashed) lines. As an external radiation field, only the 
present-day (2.73K) CMB is included although its effect can be neglected except for the lowermost 
temperature regime in the [Z/H]=0 case. We refer to this set of models as the hducial cases 
hereafter. The dotted lines indicate those of constant Jeans mass. The dashed line labeled “rj = 1” 
shows the location where the central part of a cloud becomes optically thick to continuum radiation 
(eq. 20; see discussion later). Before this condition is met, i.e. to the left of the line, clouds are 
still optically thin to continuum. 

The evolution of clouds with metallicity [Z/H] = —(4—3) at nn ~ 10^“®cm“^ is mostly affected 
by the improvements of the model, namely by the inclusion of HD cooling and the modification 
of the collapse timescale. Conversely, both at lower and higher metallicity values, the thermal 
evolution is hardly altered. 



In Figure 2, we illustrate separately the contribution by each processes to cooling/ heating 
rates in the fiducial cases. For the same cases, we also show the evolution of H 2 and HD fractions 
in Figures 3 and 4, respectively. In terms of major coolants, the thermal evolution can be classified 
into the following three metallicity ranges: (i) quasi-primordial clouds {[Z/H] < —6) (ii) metal- 
deficient clouds (—5 < [Z/H] < —3) (iii) metal-enriched clouds (—2 < [Z/H]). For each of these 
ranges, we now describe the evolution presented in Figures 1-4. 

1. Quasi-Primordial Clouds {[Z/H] < —6) 

The presence of metals at metallicity levels as low as [Z/H] < —6 does not significantly 
alter the thermal evolution in any density range. The evolution of both temperature and 
chemical species follows closely those of the metal-free case. Molecular hydrogen is always 
an important cooling agent. HD hardly affects the overall evolution despite contributing the 
cooling as much as H 2 at ~ 10^cm“^ (see also Bromm et al. 2002). 

To clarify the reasons for this trend, we first describe the HD formation process. The abun¬ 
dance of HD is determined by the equilibrium between the formation reaction (D4 reaction 
in Table 1) and its inverse dissociation reaction (D6), 

D4, D6 : D+-F H 2 ^ H+-F HD. (12) 

The deuterium ionization degree is set by the equilibrium between D1 and D2 reactions, 

D1,D2 ; D-FH+^ D+-hH. (13) 


The HD to H 2 ratio is then 

n(HD) _ /cD 4 n(D+) _ kDikDin{T)) 421K n(D) 

n{n2) ~ kD^n{H+) ~ kD^kD2n{n) ~ T ^n(H)’ ^ ^ 

where we have used the reaction coefficients given in Table 1. From equation (14), we see 
that the HD is formed abruptly if the temperature is reduced significantly with respect to 
the value corresponding to the energy difference between H 2 and HD (421K). Specifically this 
occurs when the temperature falls below about 150K. 

In the case of quasi-primordial clouds, HD becomes an important coolant at ~ 10^cm“^, after 
H 2 has already reached LTE. Almost at the same time HD reaches LTE, cooling becomes 
inefficient, and the temperature increases with increasing density. As a result of temperature 
increase, HD is dissociated immediately (see Figure 4). This is the reason why HD does not 
affect the thermal evolution in the lowest metallicity cases. 

H 2 is formed via the H~ channel (combination of reactions H2 and H3) in the low density 
regime < 10®cm“^. Because of recombinations of electrons, the catalysts of this reaction, 
the H 2 fraction saturates. This results in the plateau in H 2 fraction /h 2 ~ 10“^ in Figure 
3. Hydrogen becomes fully molecular only after the three-body reaction becomes important 
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at densities higher than the threshold value nn 10®cm ® (Palla et al. 1983). The rapid 

increase in H 2 fraction for h-h ^ 10®cm“® is owing to the three-body reaction (Figure 3). 

During the active phase of three-body H 2 formation, the heating is dominated by liberation of 
the chemical binding energy (Figure 2 a,b). It exceeds the compressional heating by about an 
order of magnitude. After the hydrogen becomes almost fully molecular at lO^^cm ®,theH2 
formation heating rate drops. However, almost simultaneously the H 2 line becomes optically 
thick and the cooling rate for a give temperature drops. Then the temperature continues to 
increase gradually. In the [Z/H]=-6 case, the grain thermal emission dominates a short interval 
around 10^®“^^cm“® before dust evaporization at ~ 10^^'^cm“® (Figure 2 b). Because of 
increasing temperature, the H 2 is about to dissociate at ~ 10^^cm“®. Indeed, a small fraction 
of H 2 dissociates. At this time, H 2 collision-induced emission (CIE) becomes important 
since the high density leads to frequent collisions and thus higher emissivity. CIE and its 
inverse collision-induced absorption are emission and absorption, respectively, of photons 
from short-lived “supermolecular” species (i.e., H 2 — H 2 ,H 2 — He, or H 2 — H) formed during 
collisions between H 2 and other particles. These supermolecules can emit and absorb photons 
with much higher probability than an isolated H 2 molecule because of their temporarily 
induced electric dipoles (e.g., Lenzuni, Chernoff, & Salpeter 1991). Collision-induced emission 
dominates the cooling until it becomes optically thick to the self-absorption at lO^^cm 
During this interval, the H 2 that was partially dissociated is re-assembled in a fully molecular 
form. The alternating behavior of H 2 formation cooling and dissociation heating in Eigure 2 
(a),(b) reflects this history. 

After the cloud becomes optically thick to H 2 collision-induced absorption, radiative cooling 
rate decreases rapidly (see Figure 2 a, b). However, the temperature continues to increase 
only gradually owing to the effective cooling by H 2 dissociation. Because the temperature is 
already very close to the dissociation temperature, the dissociation starts just after the clouds 
becomes optically thick. The temperature begins to increase almost adiabatically when the 
dissociation is nearly completed at ~ 10^®cm“®. Finally, with some more compression owing 
to inertia, a hydrostatic core, or in another word, a protostar is formed at 10^^cm“® with a 
mass-scale of ~ I0 “®Mq (Omukai & Nishi 1998). This will eventually grows to a massive 
main sequence star owing to the subsequent accretion of ambient matter (Omukai & Palla 
2001, 2003; Tan & McKee 2004). 

2. Metal-deficient Clouds (—5 < [Z/H] < —3) 

For prestellar clouds in this metallicity range, the effects of metals become more visible. 
Although H 2 -line emission still dominates the cooling in the low-density regime (< 10^cm“® 
for [Z/H]=-5; < 10^cm“® for [Z/H]=-4; < 10®cm“® for [Z/H]=-3; see Figure 2 c-e), an 
enhanced H 2 fraction due to grain surface reaction (reaction H8; see Figure 3) results in a small 
deviation in temperature evolution from the metal-free case (Figure 2). Below ~ 10®cm“®, 
differences in the H 2 fractions at comparable densities between the metal-free case and other 
cases are a result of grain-surface reaction. Because of this reaction, the H 2 fraction continues 
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to increase towards higher densities and the plateau in H 2 fraction /hj ~ 10 ^ disappears. 

The next evolutionary phase is dominated by HD cooling. For [Z/H\ = —5, this effect 
is not yet apparent, but for [Z/H]=-4 and -3, HD becomes the dominant coolant in the 
range ~ 10^ — 10®cm“^ for [Z/H]=-4 and ~ 10^ — 10^cm“^ for [Z/H]=-3 (see Figure 4 c-e). 
Unlike for quasi-primordial clouds, HD becomes important because additional H 2 formed by 
grain-surface reaction enhances cooling and makes the temperature to fall below the critical 
value ~ 150K for HD formation. Once this condition is satisfied, dramatic HD formation 
and resultant HD cooling further lowers the temperature. HD formation is thus accelerated. 
Also, unlike in the quasi-primordial clouds, the HD is not destroyed at densities > 10^cm“^ 
because for [Z/H]=-4 and -3 the temperature remains low (Figure 4). 

At higher densities, the dominant cooling mechanism is hydroxyl and/or water line emission. 
As shown in Figure 6 , the major form of oxygen compounds depends on the temperature and 
thus on the metallicity. Namely, for [Z/H]=-5 and -4, O is converted to H 2 O through OH, 
while for [Z/H]=-3 a large part of O becomes O 2 . The relevant oxygen reactions are shown 
schematically in Figure 5 (a). OH is mainly formed via the reaction 

ZIO ; 0 + H^ OH + 7 . (15) 

For [Z/H]=-5, OH is immediately converted to H 2 O by the reaction 

Z12 : H 2 OH ^ H 2 O -h H (16) 

at a density r\j lO^cm For [Z/H]=-4, OH is formed at lO^cm ^ and remains in this form 
for some intervals because of the lower temperature with respect to the [Z/H]=-5 case. Some 
fraction of OH becomes O 2 via the reaction 

Z19 : O + OH^ O 2 + H. (17) 

The reaction coefficient of the H 20 -bearing reaction Z12 is kzi 2 oc r^'®exp(—1660iF/T) while 
that of the 02 -bearing reaction Z19 is kzi 9 oc r“®'^exp(—30A'/r). Then, in a low tempera¬ 
ture environment, H 2 O formation is significantly suppressed and O 2 is the dominant oxygen 
compound with respect to H 2 O. In particular, for T < 300K, the reaction Z12 is frozen and 
O is locked in O 2 without forming H 2 O. At densities ~ 10 ®cm“^, the temperature rises sud¬ 
denly owing to heat injection by the three-body H 2 formation. Consequently, OH and also 
O 2 become H 2 O via the reactions Z20 and Z12, respectively. For [Z/H]=-3, the temperature 
remains low (< 200K) and the majority of O becomes O 2 and remains in this form without 
changing to H 2 O. 

During this OH/H 2 O cooling phase, H 2 formation heating exceeds the compressional heating 
(Figure 2 c-e). As a consequence, the temperature increases at < 10®cm“^ for [Z/H]=- 
5, (< 10®cm“^ for [Z/H]=-4, and < 10^cm“^ for [Z/H]=-3, respectively). For [Z/H]=-5, 
the temperature is sufficiently high to hide this effect which instead becomes apparent for 
[Z/H]=-4. When [Z/H]=-4, heat injection due to the abrupt onset of the three-body reaction 
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causes a dramatic temperature rise at ~ 10®cm“^ since the temperature in the preceding 
evolution is low. On the other hand, for [Z/H]=-3, heat injection is more gradual since the 
majority of H 2 is formed by the more gradual grain-surface reaction. 

When H 2 formation dominates the heating mechanism, the major cooling agent is represented 
by dust thermal emission (see Figure 4c-e). Dust and gas exchange energy by collisions. At 
high density, frequent collisions lead to a higher cooling rate (A oc n per unit mass). In 
spite of this cooling, the temperature continues to increase as long as H 2 formation heating 
is important. 

When H 2 formation is completed, the heating rate drops, and the temperature begins to fall 
rapidly due to dust cooling. Dust cooling remains dominant until the clouds become optically 
thick to dust self-absorption, at densities in the range 10^^“^^cm“^ for —5 < [Z/H] < —3. 
After that, the cloud contracts adiabatically and temporary hydrostatic cores (first cores) 
might be formed as for present-day star formation (Larson 1969; see below). When the 
temperature reaches 2000K, H 2 dissociation occurs. The evolution after this point is exactly 
the same as that of the quasi-primordial clouds. 

3. Metal-enriched Clouds (—2 < [Z/H]) 

When metal-enrichment proceeds up to [Z/H] ~ —2, the thermal evolution of prestellar 
clouds becomes qualitatively the same as the present-day one. Series of coolants become 
dominant one after the other in the same order. Only the characteristic density of each 
phase is shifted towards lower densities at higher metallicities. Note the similar shapes in 
the temperature-evolution curves in this metallicity range in Figure 1. The rapid initial 
temperature drop is induced by fine-structure line emission of C'*', C and O. Unlike in the 
metal-deficient clouds, HD cooling does not affect the thermal evolution since C, O cooling 
are dominant. The temperature increase (at densities 10^“®cm“^, 10^“^'^cm“^ for [Z/H]=-2 
and -1, respectively) and stagnation (at 10*^“^cm“^ for [Z/H]=0) are due to H 2 formation 
heating although the grain-surface reaction deposits only a few % of the binding energy (0.2 
eV of the total 4.48 eV in the lowest density regime); the rest of the energy is used for exciting 
H 2 ro-vibrational lines or deposited on the grain and eventually radiated away in both cases 
(Hollenbach & McKee 1979). 

In the meantime, atomic carbon is converted to carbon monoxyde (CO). Major carbon-related 
reactions, including CO formation, are shown schematically in Figure 5 (b). CO cooling is 
dominant in the density range 10^'® — lO^'^cm"^ (10^ — 10®cm“^) for [Z/H]=-2 ([Z/H]=-l) 
although for [Z/H]=0 it is significant only for a short time interval around 10^cm“^ (Figure 2 
f-h). In some other studies, CO can be the dominant coolant for a wider range of densities even 
for [Z/H]=0 (e.g., Koyama & Inutsuka 2000). Probably this is due to different assumptions 
on the external radiation field and on the effects of cosmic rays. At higher densities, infrared 
thermal emission by dust grains becomes important. The gas temperature decreases until 
it reaches the dust temperature, and then thermally couples to dust grains. This occurs at 
densities ~ 10®cm“^ for [Z/H]=-2, ~ 10®'^cm“^ for [Z/H]=-l, and ~ 10®cm“^ for [Z/H]=0. 
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After that, the thermal evolution is controlled by dust cooling and the temperature increases 
gradually with increasing density. When the central parts of the clouds become optically 
thick to dust self-absorption, the temperature begins to increase adiabatically. The evolution 
thereafter is the same as in the previous case. 


After the clouds become optically thick to continuum, the temperature increases adiabatically. 
As seen in Figure 1, all the paths in the adiabatic contraction converge approximately to the dashed 
line in spite of different values of temperature and density when the clouds become optically thick. 
The reason for this convergence has already been explained in Paper I, but we repeat it here briefly 
for later use. First, let us consider the temperature evolution paths just before the clouds become 
optically thick. Since cooling is dominated by continuum emission (i.e., the H 2 collision-induced 
emission or the dust thermal emission depending on metallicity) and the heating is dominated by 
compressional heating, the temperature is determined by the balance between them: 

AkzctT* = (18) 

^col 

where kz is the continuum opacity for gas with metallicities Z and we use the pressure-free free-fall 
timescale tcoi = ^coi.o = \/2>'al‘i2Gp for this analytic estimate. Note that the temperature evolution, 
T = r(nH), as given by equation (18), depends on metallicity through the opacity nz- While the 
density increases along this path, the central region, whose length scale is ~ Aj, becomes optically 
thick at the position in the nn — T plane where the condition 


TJ = KzXj = 1 

is satisfied. Eliminating kz from equations (18) and (19), we obtain the relation 


T = 






H 


(19) 


( 20 ) 


This is shown by the dashed line denoted as “rj = 1” in Figure 1. So far, the meaning of this line 
has been that of the line connecting all the points where the central parts of the clouds become 
optically thick for different metallicities. However, this line has another significance. Note that the 
exponent of nn (=2/5) in this relation is accidentally equal to 7ad — 1 — 0.4 for molecular hydrogen. 
Then, this line coincides with an adiabatic path T oc Since this path passes through the 

points where the clouds become optically thick, this is indeed the adiabatic evolutionary line in 
optically thick phase for any Z. This is the reason for the convergence of the adiabatic paths. 
In other words, when the clouds becomes optically thick, their specific entropy has the universal 
value. The evolution in the optically thick phase follows just the constant entropy path with this 
universal adiabat, if entropy in chemical binding energy is taken into consideration. 


3.1.2. Effects of Cosmic Microwave Background Radiation 

Next we discuss the effect of an external radiation field on prestellar thermal evolution. In this 
paper we consider only the CMB in high-redshift universe. Figure 7 shows the cases irradiated by 
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the black-body radiation with Trad = 54.6iir corresponding to 1 + z = 20, the typical epoch for the 
first star formation in the universe. The cases without strong CMB (fiducial cases shown in Figure 
1) are indicated by dotted lines for comparison. 

As expected, the temperature cannot drop below Trad- On the other hand, the behavior at 
higher temperatures is hardly affected by the CMB radiation. An exception to this is the case of 
[Z/H]=-3 around ~ 10®cm“^, where the temperature jumps up owing to heat injection by three- 
body H 2 formation in the strong CMB case. In this case, as a result of suppression of H 2 formation 
on grains by higher dust temperature, three-body reaction becomes more important than in the 
fiducial case. 

Another difference is that the temperature in the adiabatic contraction phase does not coincide 
with the universal adiabat. The reason is clear from the discussion above, where we have assumed 
that heating is dominated by compressional heating when the clouds becomes optically thick. For 
higher metallicities, this condition is violated because of CMB heating of dust grains. Despite the 
non-convergence in the adiabatic phase, after this phase these evolutionary paths converge to that 
of the zero-metallicity case. Indeed, the evaporation of dust grains suddenly decreases the optical 
depth in the central parts of the clouds. After that, the clouds evolve by H 2 CIE cooling, exactly 
in the same as in the zero-metallicity case. This is partly an artifact of our one-zone treatment 
because in the envelope we expect dust grains to survive and to remain optically thick even though 
the central parts has become optically thin again. In reality, photons emitted in the central region 
are absorbed in the envelope that is heated and that begins to back-warm the central region. If this 
effect is important, our treatment breaks down. To explore this, a radiation hydrodynamic study 
in necessary is required. This is beyond the scope of this paper, and we postpone it for a future 
study. 


3.1.3. Role of Metals in Gas and Dust Phases: A Minimum Model for Metallieity Effects 

Metals can be either in gas phase or dust grains. These two forms of metals influence the 
thermal evolution in different ways. Grains are of importance to (i) cooling by dust grains, and 
(ii) H 2 formation on grain surface, while gas-phase metals act via (iii) cooling by atomic/molecular 
metal lines. To see these effects separately, we show the same cases as the hducial ones (Figure 1) 
but without gas-phase metals in Figure 8, and without dust grains in Figure 9. Other parameters 
are the same as in the hducial cases, whose temperature evolution are shown by dotted lines for 
comparison. Although we still use the same denomination [Z/H] for easier comparison with the 
hducial cases, the amount of total metals are smaller than that in the hducial cases with the same 
[Z/H] because one form of metals is missing. 

Figure 8 shows that with dust grains only, the thermal evolution for gas metallicities [Z/H] 
< —3 is approximately reproduced. For higher gas metallicities ([Z/H] > —2), the behavior in the 
low density regime is strongly altered because of the lack of the hne-structure-line cooling. In the 



- 14 - 


metallicity range [Z/H]~ -(5-3), cooling in the intermediate density regime is dominated by OH 
and H 2 O, which have not been included in this Figure. Despite this, there is only a rather small 
deviation in temperature since cooling by H 2 and HD have comparable effects. Of the two effects of 
dust grains, dust cooling becomes important at lower metallicity than the effect on H 2 formation. 

Figure 9 shows the predicted thermal evolution when no dust is present. With gas-phase metals 
only, the temperature behavior is completely altered. From the comparison between Figures 8 and 
9, we can conclude that at low metallicities cooling by dust grains is dominant and the thermal 
evolution can be reproduced even without gas-phase metals. On the other hand, in metal-enriched 
cases, cooling by fine-structure lines are important in the low density regime. 

Considering only the dust and the fine-structure cooling, we can construct a minimum model 
of metallicity effects. In addition to the metal-free model which includes H 2 and HD cooling, 
H 2 collision-induced emission/absorption, and Hl- 8 , D-1-6 reactions, we include fine-structure line 
cooling and recombination of C'*' (Z28). The result is shown in Figure 10. Indeed, this model 
provides a good approximation of the thermal evolution. Although there are still small deviations 
owing to the lack of OH and H 2 O cooling, this model significantly simplifies the implementation 
of metallicity-dependent effect, since it does not solve the complex metal chemistry nor the level 
populations of C and O molecules. As the fine-structure lines are two (CH) or three (Cl, 01) level 
systems, we can easily include them. 


3.2. Reduced Chemical Model for Low-metallicity Clouds 

Although the minimum model above succeeds in roughly reproducing the thermal evolution 
of our full model, there are still some deviations. As a further step, here we develop a model that 
reproduces the temperature evolution as well as the abundances of the dominant chemical species. 

In the full model, we have incorporated about 500 reactions (478 reactions in paper I and 18 
deuterium reactions) between 50 species (H, H 2 , e“, H'*', H^, H)|“, H“, He, He^, He^"'', HeH"'', D, 
D+, D-, HD, HD+, C, C 2 , CH, CH 2 , CH 3 , CH 4 , C+, C^, CH+, CH+, CH+, CH+, CH+, O, O 2 , 
OH, CO, H 2 O, HCO, O 2 H, CO 2 , H 2 CO, H 2 O 2 , 0+, O^, OH+, CO+, H 2 O+, HCO+, O 2 H+, H 3 O+, 
H 2 CO'’', HCO^ and H 3 CO^). Apparently, not all of these reactions are of equal importance. In 
many cases, we are only interested in reproducing the thermal history accurately and we are not 
interested in reproducing the evolution of all chemical species correctly. To this end, we need to 
properly account for the abundances of the dominant coolants. We select a reaction network which 
reproduces the abundance of dominant coolants, as well as dominant H, C and O compounds and 
the electron fraction in any evolutionary stage before the prestellar clouds becomes optically thick. 
This model includes 23 species H, H+, H", H 2 , e, D, D+, HD, O, OH, H 2 O, O 2 , 0+, OH+, H 2 O+, 
H 30 ''', O^, C"'', C, CH, CH 2 , CH 3 , CH 4 , CO, CO"'', CO 2 and 54 reactions, which are given in Table 
1. Reactions related to C and O are shown schematically in Figure 5. 

In Figure 11, we present the temperature evolution predicted by the model with full network 
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and that by the reduced one. The agreement is so good that the curves appear to be degenerate in 
the plot. The corresponding evolution of the dominant chemical species are also shown in Figure 12, 
for comparison. Although there are still some deviations when the fractional abundances are low, 
they are not important coolants at those times. When some species are either important coolants 
or represent the dominant compound of a given element, their abundances are well reproduced by 
the reduced model. 


4. Summary and Discussion 

In this paper, we have studied the thermal and chemical evolution of star-forming clouds with 
different metallicities. We assume that the collapse rate of the clouds is approximately at the free- 
fall rate and the characteristic radius of collapsing clouds is given by the Jeans length. The model 
developed in Omukai (2000) has been updated under some aspects, most notably by inclusion of 
deuterium chemistry. The temperature evolution obtained from this model has been presented in 
Figure 1. The HD cooling modifies evolution of the clouds in the metallicity range Z 10-(3-5)Zq, 
this species being the dominant coolant in the density range 10^ ®cm 

While we have solved the thermo-chemical microphysics in great detail, we have adopted a 
simplified dynamical evolution. Our full thermal modeling is numerically too costly to be coupled 
to detailed hydrodynamical models. To favor such experiments, though, a reduced model retaining 
the essential features of the full model is required. As a simplest model of this kind, we have shown 
that a model including only fine-structure line cooling and dust process, i.e. dust-cooling and H 2 
formation on dust surface, in addition to the metal-free gas physics, reproduces the temperature 
evolution very well. A better (but more time-consuming) model has also been implemented with 
reduced chemistry (with about 20 species and 50 reactions) which can reproduce the thermal 
evolution as well as the chemical evolution of important species. 

In Appendix, by applying the linear theory of cloud cores deformation to the above temperature 
evolution, we have discussed fragmentation history of the star-forming clouds. Given the initial 
distribution of elongation, the mass spectrum of fragments can straightforwardly be derived from 
this theory. As examples, we have shown the mass spectrum for two different initial distributions 
of core elongations: (i) a uniform initial distribution, and (ii) a Gaussian distribution peaking 
around half of the critical value of fragmentation with variance of 0.1. In both cases, we have 
found that while in the metal-free case fragmentation is possible only at the end of the efficient 
molecular cooling phase, fragmentation at high density < 10^*^cm“3 ig enabled even with a slight 
—6) enrichment level via dust grains, the probability of this fragmentation increases 
drastically around [.Z’/Ff] ~ —5. Since the fragmentation mass-scale is determined by the thermal 
Jeans scale at the fragmentation epoch, this results in sub-solar mass fragments even for slightly 
metal-enriched cases. 

According to observations of nearby star-forming regions, stellar initial mass function (IMF) 
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does not appear to depend strongly on metallicity (e.g., Kroupa 2002). However, the metallicity 
range studied is still very limited (a factor of a few). We predict that if systems with metallicities 
around [Z/H] ~ —4 could be sampled, their IMF should resemble that from the present results. In 
the quasi-primordial case [Z/H] < —6, very massive fragments are predominantly produced. Even 
in this case if extremely strong UV field is present photodissociation might lead to low-mass star 
formation by atomic cooling (Omukai & Yoshii 2003), although Bromm &: Loeb (2003b) argue an 
alternative possibility that a supermassive object is the likely outcome for the atomic-cooling cloud. 
For modestly metal-enriched {[Z/H] < —1) clouds, fragmentation can occur at two different stages 
during the collapse, i.e. both in the molecular-cooling phase at lower density and in the dust-cooling 
phase at higher density. The former fragmentation produces predominantly massive stars, although 
the typical mass-scale decreases with increasing metallicity; the latter produces sub-solar mass stars. 
The relative importance of these two fragmentation modes changes with metallicity. Dust-induced 
fragmentation is of increasing importance for clouds with higher metallicity. Reflecting the two 
modes of fragmentation, the IMF often shows two peaks in low-metallicity star-forming clouds. 
As metal enrichment increases to the level of [Z/H] ~ —1, these two peaks eventually merge as 
the line- and dust-cooling phases become indistinguishable in temperature evolution. This might 
explain why the present-day IMF is single-peaked. 

In this work, we have assumed the same dust-size distribution as the standard MRN distri¬ 
bution for local ISM (Mathis, Rumpl, & Nordsieck 1977). The origin of dust grains in the local 
ISM is largely attributed to asymptotic giant branch (AGB) stars with smaller contribution from 
Type II SNe. In the early universe dust mainly comes from Type II (Todini & Ferrara 2001) or 
pair-instability SNe (Nozawa et al. 2003; Schneider, Ferrara & Salvaterra 2004), since there could 
not be sufficient time for star to evolve to the AGB phase. The average SN grain radius is smaller 
than for AGB stars, and thus the total area per unit mass of dust is larger than the standard 
MRN distribution used in this work. The large grain area enhances both dust-gas thermal coupling 
and the H 2 formation rate on grain surface. Therefore, even small dust amounts can produce pro¬ 
nounced effects on the chemistry. Also the elemental abundance ratios in ejecta from the first SNe 
are different from those in local ISM, which we have used in the analysis in this paper. Elements 
important for the cooling are carbon and oxygen. In comparison with the local ISM, carbon is 
significantly deficient (< 1/100) in the first SNe ejecta while the amount of oxygen is similar in 
both environments. Among these two elements, in low metallicity environment atomic and molec¬ 
ular forms of oxygen are more important coolants than those of carbon (§3). As gas-phase oxygen 
abundance is similar, difference in gas-phase abundance ratio does not alter thermal evolution of 
low metallicity clouds significantly. We plan to study effects of different grain size distribution and 
gas phase elemental abundances in a forthcoming paper. 

Recently, Bromm &: Loeb (2003a, BL03) have also discussed the critical metallicity separating 
low- and high-mass star formation modes. In their analysis, they included only fine-structure-line 
cooling of G"*" and O, and derived the critical metallicity by estimating the amount of metals above 
which the fine-structure-line cooling becomes more efficient than the compressional heating (at 
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density ~ 10^cm“^) where the temperature drops to a minimum (~ 200K) and fragmentation is 
expected to occur. For cases in which either [CII] or [01] is present, they obtained critical metallicity 
values of [C/H]cr\t = —3.5 and [0/H]cnt = —3.05, respectively, by comparing the cooling rates at 
10^cm“^ and 200K. The comparison between their results and those derived in this work is shown 
in Figure 13, where we plot the temperature evolution calculated by including only (a) [CII] or (b) 
[01] cooling in addition to the metal-free thermo-chemical processes. In the case (a) in which only 
[CII] cooling is included, the position of temperature minimum is consistent with BL03 analysis. 
For the (b) [01] case, however, despite BL03 claim that [0/i7]crit = —3, even with [0/H] = —4 
the position of the temperature minimum is already shifted to nn ~ 10®cm“^. This is because [01] 
cooling becomes important only at densities higher than the value 10^cm“^ adopted by BL03. Note 
that the temperature at nn ~ 10^cm“^ is hardly affected. In our model, in which dust grains as 
well as gas-phase metals are considered, at [Z/H] = —4, the temperature becomes much lower due 
to the combined effects of grain-surface H 2 formation and the HD cooling. Thus, in the absence of 
dust cooling, the fragmentation mass-scale remains relatively high, > IOMq, and sub-solar mass 
stars are not formed. Hence we conclude that the dust-induced fragmentation leads to low-mass 
star formation for metallicity values well below those predicted by BL03. 


We thank the referee, Volker Bromm, for comments that improved the presentation of the 
paper. This study is supported in part by the Grants-in-Aid by the Ministry of Education, Science 
and Culture of Japan (16204012:KO,14740129:TT). 


A. Fragmentation Mass-scale 

In this section, we present a semi-analytic description of the hydrodynamical evolution for 
prestellar clouds whose thermal evolution has been described above. In particular, we discuss the 
fragmentation history allowing to estimate the characteristic mass-scale of the fragments. 

Following the results of numerical simulations of the collapse of metal-free prestellar clouds 
(e.g., Bromm et al. 1999; 2002), we expect that these clouds experience fragmentation at the first 
temperature minimum, i.e. when the effective adiabatic exponent dlogp/dlogp = 7 exceeds 1 for 
the first time after the rapid cooling phase. For metal-free clouds, this happens at number densities 
nn ~ 10^cm“^. After this fragmentation, the clouds start the run-away collapse, which can be 
described by the Larson-Penston-type self-similar solution, which was originally found by Larson 
(1969) and Pension (1969) for the isothermal case, 7 = 1, and was extended for other values of 7 < 
4/3 by Yahil (1983). We assume that collapsing prestellar clouds follow the Larson-Penston-type 
self-similar solution with the effective adiabatic exponent 7 as derived from the thermal evolution 
discussed above. 

In particular, we discuss whether run-away-collapsing clouds fragment again or not. Linear 
theory on non-spherical deformation of the self-similar solution has been developed by Hanawa & 
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Matsumoto (2000) and Lai (2000). According to these authors, the amplitude of a perturbation 
grows in a power-low of density; if we define the elongation as £ = {b — a)/a where a, b are minor 
and major axis lengths of the cores, then £ cc with n being a constant for a given 7 . The 
power of growth n is presented in Figure 14 as a function of 7 . As Hanawa & Matsumoto (2000) 
and Lai (2000) considered only values of the adiabatic exponent in the range 0.9 < 7 < 1.25; 
outside this range of values, we extrapolate n in the following way. For 7 > 1.25, the elongation 
decreases rapidly with increasing density and we linearly extrapolate the power of growth n from 
lower values of 7 . When 7 is reduced below 0.9, the solution approaches the homogeneous density 
distribution with the density contrast between the center and the envelope decreasing. Eventually, 
at 7 = 5/6, the solution becomes the uniform-density one, for which the growth rate of elongation 
is n = 0.354. Here, we use the linear interpolation for 5/6 < 7 < 0.9 and for 7 < 5/6 we use that 
for the uniform-density solution. Our treatment of assuming the self-similar solution breaks down 
when 7 > 4/3, where no self-similar solution exists. In reality, when 7 exceeds 4/3 at some density, 
a hydrostatic core forms at the center and it contracts owing to the accretion of the envelope after 
that. This phase is treated by the same prescription as in §2. 

When the elongation grows to a certain level of non-linearity, fragmentation is expected to 
occur since filamentary clouds are gravitationally unstable (e.g. Larson 1985). Given the lack 
of systematic investigation of the non-linear evolution, we use a simple prescription, neglecting 
complicated effects of, e.g., rotation etc. Initially, just after the first fragmentation, prestellar 
clouds are assumed to be elongated with amplitude Ti(< 1). At this point, a functional form 
for the probability distribution for £[ is assumed. As the clouds collapse, the elongations grow (or 
decay) following the linear theory. When the elongation exceeds a threshold value Tnl ~ 1; the cores 
becomes filamentary and those filaments fragment into two cores. The mass of fragments are given 
by the thermal Jeans mass at the fragmentation. Although the cores might grow further by accretion 
of the surrounding gas after fragmentation, this effect is not included. The elongation of newly 
produced fragments obeys the same probability distribution as that of the initial cores. Evolution 
of those fragments thereafter is treated in the same way as the initial cores. By repeating the above 
procedure, we can obtain a mass distribution of fragments. Since the probability distribution of 
elongation is unknown, we make two different assumptions: (i) a uniform initial distribution of 
elongations and (ii) a distribution peaking around half of the critical value of fragmentation with 
variance a = 0.1. For each case, we have carried out Monte-Carlo simulations with sample number 
8192. By reducing the sample number to half this value, we have confirmed that this is sufficient 
for convergence. 

In Figure 15, we show the evolution of elongation for the fiducial cases. If the elongation 
exceeds the threshold Tnl at some point, the clouds are assumed to fragment and the elongation 
is reset to an initial value. This effect is not included in this Figure. Suppose this occurs at 
density nH,frag and a new fragment is assigned with a value of elongation Tiyrag- At higher density 
> nn.frag, the ratio T/Ti,frag is given by (T/Ti)|„jj/(T/Ti)|„jj using £/£i shown in the Figure. 

In Figure 16, the mass distribution of fragments is shown by histogram for each metallicity 
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case. The two panels correspond to the different initial distributions of elongation. The results 
in the two cases are qualitatively very similar. As our model is still rather crude, we should not 
take these quantitative differences seriously and should focus only on qualitative agreement. In the 
metal-free case, only massive (~ 10 ^Mq) fragments are produced without fragmenting in the run¬ 
away collapse phase. This can be understood from Figure 14. The maximum value of elongation 
is attained around nn ~ 10^®cm“^, but it is only as large as the initial value Si and then less than 
Tnl- In contrast, for metallicities as low as [Z/H]=-6, low-mass fragmentation is possible although 
the fraction in such objects is still small. The reason is that since the maximum elongation exceeds 
the initial value, the core can fragment provided that the initial elongation is close to the critical 
value Tnl- Since highly elongated cores are more numerous for the uniform Si distribution, the 
fraction of low-mass fragments is higher in this case. When metallicity reaches [Z/H]=-5, the 
low-mass stars (~ O.IMq) outnumber very massive objects (~ W^Mq) because the maximum 
elongation (around lO^^cm”^) is sufficiently higher than the initial value and then the elongation 
can grow sufficiently around 10^^cm“^. In other words, the temperature minimum by the dust¬ 
cooling becomes sufficiently deep and fragmentation at this point occurs with high probability. In 
the [Z/H]=-4 case, as a result of the sudden heat-up around 10®cm“^, elongation of cores decreases 
rapidly (see Figure 15). Without seed elongation, the fragmentation does not occur at higher 
density even with dust cooling. Therefore, fragmentation occurs only at the line-cooling phase 
and produces only rather massive ~ IOMq fragments. Heat injection being gradual, the situation 
in the case of [Z/H]=-3 is rather similar to that of [Z/H]=-5. Because fragmentation occurs at 
10^® — lO^^crn”^ with high probability (see Figure 15), only a single peak of low-mass objects 
appears in the mass distribution. For [Z/H]=-(2-l), two peaks appear in the mass distribution 
corresponding to two minima of the temperature evolution (see Figure 1): one is around 10 ^Mq, 
resulting from line-induced fragmentation around 10^“^cm“^, and the other around IMq resulting 
from the dust cooling around 10®“'^cm“^. With [Z/H]=0, only a single peak appears IMq 
since the two temperature minima almost merge into a single one (Figure 1). In metal-enriched 
cases, because of the lower temperature, the thermal evolution is more vulnerable to external 
heating effects than in lower metallicity cases. External heating by stellar/active-galaxy radiation 
or cosmic rays is quite uncertain in primeval galaxies and it is not included in our calculation. 
Because of this, results concerning the high-metallicity cases are to be taken with some care. With 
this caveat in mind, resemblance between our fragment mass distribution and the observed IMF is 
encouraging. 

Our simple analysis did not take into account interaction between fragments and accretion after 
fragmentation. However, it is known that the IMF in present-day star formation is largely shaped 
by the fact that stars form in clusters out of gas clouds that are pervaded by supersonic turbulence, 
through the complex dynamics of many interacting and accreting protostars (e.g. Bate, Bonnel, 
& Bromm 2003). Such a clustered nature is also vital to predict the shapes of prestellar clouds 
themselves due to the gravitational interactions in the cluster. Besides, fragmentation occurs not 
only in dense molecular cores, but also in massive circumstellar disks after formation of a central 
star (Bate et al. 2003). Because of those processes not included in our analysis, fragment mass 
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function obtained here should be regarded as showing only typical mass-scale rather than detailed 
mass distribntion. In order to stndy the core mass fnnction more quantitatively in future, we need 
to consider all those cluster effects above. 
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Table 1. Reduced Chemical Reactions for Low-Metallicity Prestellar Clouds 


Number Reaction Rate Coefficient ^ Reference 

HI H+ + e^H + 7 fcni = exp[-28.6130338 - 0.72411256(ln(T(eV))) A97 

-2.02604473 x 10-2(ln(T(eV)))2 
-2.38086188 x 10-3(ln(T(eV)))3 
-3.21260521 X 10-^(ln(T(eV)))^ 

-1.42150291 X 10-5(ln(T(eV)))5 
+4.98910892 x 10-6(ln(T(eV)))® 

+5.75561414 x 10-^(ln(T(eV)))^ 

-1.85676704 x 10-8(ln(r(eV)))8 
-3.07113524 X W-^{ln{T{eY)) f] 


H2 

H + e^ 

1 

+ 

kH2 = 1.4 X 10-^^T^-^‘^^ex.p{-T/1.62 x 10^) 

GP98 

H3 

H- + H 

^ H 2 + e 

kH3 = 4.0 X (T > 300) 

kH3 = 1.5 X 10-9 (T < 300) 

GP98 

H4 

H 2 + H- 

3H 

kHA = kl^'^kl 





/cl = 1.12 X 10-i9exp(-7.035 x lOV^) 

SK87 


kn = 6.5 X 10-'^T-^/2 

exp(-5.2 X 10^/r)[l - exp(-6000/r)] 
a = (1 + n/ncr)~^ 

logio(ncr) = 4.0 - 0.4161ogio(r/10^) 

-0.327(logio(T/104))2 


H5 

3H ^ H 2 + H 

kH5 

= 5.5 X io-29r-i 

PSS83 

H6 

2H + H 2 —> 2 H 2 

km 

= km/'^ 

PSS83 

H7 

2 H 2 ^ 2H + H 2 

km 

^high +OW 

SK87 


how = 1-18 X 10 ^°exp(— 6.95 x 10^/T) PSS83 (fehigh) 

khigh = 8.125 X l0-8r-V2 

exp(-5.2 X 10^/T)[1 - exp(-6.0 x lOV?')] 
a = (1 + n/ncr)~^ 

logio(ncr) = 4.845 - 1.31ogio(R/10^) 

+1.62(logio(r/10^))' 

H8 2H + grain ^ H 2 kns = 6.0 x 10 -^^(T/ 300 K )^/2 XH85 

x[l + 4.0 X 10-2(T + rgr)^/2 

+2.0 X 10"3r + 8.0 X 10“®T2]-i 

! Z\oQ^\ 

/a = [1 + exp(7.5 X 102(1/75 - l/T^h)]-^ 
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Table 1—Continued 


Number Reaction Rate Coefficient ^ Reference 


D1 

D + H+ ^ 

D+ + H 

koi = 3.7 X 10-i°r°-28exp(-43/T) 

NU02 

D2 

D+ + H ^ 

D + H+ 

kj^2 = 3.7 X lo-iOr°-28 

NU02 

D3 

D + H 2 ^ 

D + HD 

A:z 53 = 9.0 X 10-i^exp(-3876/r) 

NU02 

D4 

D+ + H 2 - 

H+ + HD 

koA = 2.1 X 10-9 

NU02 

D5 

HD + H ^ 

►H 2 + D 

kD5 = 3.2 X 10-i^exp(-3624/r) 

NU02 

D 6 

HD + H+ - 

H 2 + D+ 

koQ = 1.0 X 10-9exp(—464/T) 

NU02 

Z1 

0 + + H^ 

H+ + 0 

kzi = 6.80 X 10-19 

M97 

Z2 

H+ + 0 ^ 

0+ + H 

kz 2 = 7.00 X 10-i9exp(-232/T) 

M97 

Z3 

0 + + H 2 - 

0H+ + H 

kz'i = 1.70 X 10-9 

M97 

Z4 

0H+ + H 2 

^ H 2 O+ + H 

kzA = 1-01 X 10-9 

M97 

Z5 

H 2 O+ + H 

2 ^ H 3 O+ + H 

kz^ = 8.30 X 10-19 

M97 

Z 6 

H 2 O+ + e 

^ O + H 2 

kzQ = 2.00 X 10-^(T/300)-9-50 

M97 

Z7 

H 2 O+ + e 

^ OH + H 

kz7 = 1.60 X 10-'^(T/300)-9-50 

M97 

Z 8 

H 3 O+ + e 

^ OH + 2H 

kzs = 6.50 X 10-'^(T/300)-9-50 

M97 

Z9 

H 3 O+ + e 

^H 20 + H 

kz9 = 3.50 X 10-^(T/300)-9-50 

M97 

ZIO 

0 + H ^ OH + 7 

kzw = 9.90 X 10-19(T/300)-0'38 

M97 

Zll 

O + H 2 ^ 

OH + H 

kzii = 3.43 X 10-i9(T/300)2-67exp(-3160/T) 

M97 

Z12 

H 2 + OH - 

H 2 O + H 

kzi2 = 1.55 X 10-i9(T/300)i-99exp(-1660/T) 

M97 

Z13 

20H ^ H 2 O + 0 

kzi 3 = 1.65 X 10-i9(T/300)i-i^exp(-50/T) 

M97 

Z14 

H+ + OH 

0H+ + H 

kziA = 2.10 X 10-9 

M97 

Z15 

H+ + H 2 O 

^ H 2 O+ + H 

kzi5 = 6.90 X 10-9 

M97 

Z16 

H + OH^ 

►H 2 + O 

kzi& = 7.00 X 10-i^(T/300)2-89exp(-1950/T) 

M97 

Z17 

H + H 2 O - 

OH + H 2 

kzn = 6.83 x 10-i9(T/300)i-99exp(-9720/T) 

M97 

Z18 

20 ^ O 2 + 7 

kzis = 4.90 X 10-20(T/300)1-®® 

M97 

Z19 

O + OH- 

^ O 2 + H 

kzi 9 = 4.34 X 10-11 (T/300)-0'i’Oexp(-30/T) 

M97 

Z20 

H + O 2 ^ 

OH + 0 

kz 2 o = 3.30 X 10-i°exp(-8460/T) 

M97 

Z21 

H+ + O 2 - 

0^ + H 

A:z21 = 2.00 X 10-9 

M97 

Z22 

O 2 + e ^ 

20 

kz 22 = 1.95 X 10-i'(T/300)-°'l’° 

M97 

Z23 

0 + CH^ 

> CO + H 

kz2‘i = 6.60 X 10-11 

M97 

Z24 

C + OH^ 

> CO + H 

kz2A = 1.10 X 10-iO(T/300)9-50 

M97 

Z25 

C + O 2 ^ 

CO + 0 

kz2b = 3.30 X 10-11 

M97 

Z26 

C+ + O 2 - 

0+ + CO 

kz2& = 6.20 X 10-19 

M97 

Z27 

OH + CO 

CO 2 + H 

kz27 = 1.00 X 10-19 

M97 
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Figures 



0 5 10 15 20 


number density log n^^ (cm“^) 


Fig. 1.— Temperature evolution of prestellar clouds with different metallicities. Those with metal- 
licities [Z/H]=-oo (Z=0), -5, -3, and -1 (-6, -4, -2, and 0) are shown by solid (dashed) lines. Only the 
present-day CMB is considered as an external radiation field. The lines for constant Jeans mass are 
indicated by thin dotted lines. The positions where the central part of the clouds becomes optically 
thick to continuum self-absorption is indicated by the thin solid line (equation 20). The intersection 
of the thin solid line with each evolutionary trajectory corresponds to the epoch when the cloud 
becomes optically thick to the continuum. To the right of this line, the clouds are optically thick 
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Table 1—Continued 


Number 

Reaction 




Rate Coefficient 

Reference i’ 

Z28 

C+ + e^ 

C + 7 

kz28 

= 4.40 

X 

10-12 ('7^/30o)-o.61 

M97 

Z29 

C+ + OH 

^ CO+ + H 

kz29 

= 7.70 

X 

10-10 

M97 

Z30 

CO+ + H 

^ H+ + CO 

kzm 

= 7.50 

X 

10-10 

M97 

Z31 

C + H ^ ' 

CH + 7 

kz3i 

= 1.00 

X 

10-17 

M97 

Z32 

C + Hs ^ 

CH + H 

kz32 

= 6.64 

X 

10-i°exp(-11700/T) 

M97 

Z33 

H + CH- 

> C + H 2 

kz33 

= 4.98 

X 

10-11 

M97 

Z34 

H2 + CH - 

CH 2 + H 

kz3A 

= 2.38 

X 

10-i°exp(-1760/T) 

M97 

Z35 

H + CH2 - 

CH + H 2 

kz35 

= 2.70 

X 

10-10 

M97 

Z36 

H2 + CH2 

^ CHs + H 

kz36 

= 5.18 

X 

10-ii(T/300)°i7exp(-6400/T) 

M97 

Z37 

H + CHs - 

CH 2 + H 2 

kz37 

= 1.00 

X 

10-i°exp(-7600/T) 

M97 

Z38 

H2 + CHs 

^ CH4 + H 

kz38 

= 6.86 

X 

10-i^(T/300)2-74exp(-4740/T) 

M97 

Z39 

H + CH4 - 

H2 + CHs 

kz39 

= 5.82 

X 

10-i3(T/300)3-00exp(-4045/T) 

M97 

Z40 

H2 + C^ 

CH2+7 

kz40 

= 1.00 

X 

10-17 

M97 


“^The temperature T is in K, except otherwise noted. 

^A97; Abel et al. (1997), GP98; Galli & Palla (1998), SK87; Shapiro & Kang (1987), M97; Millar 
et al. (1997), NU02; Nakamura & Umemura (2002), PSS83; Palla et al. (1983), TH85; Tielens Sz 
Hollenbach (1985) 
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and there is little radiative cooling. 
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Fig. 2.— Energy balance diagrams. Contribution to total cooling and heating rates (per unit 
mass) from each processes are shown as a function of number density nn for fiducial cases with 
(a)[Z/H]=-oo, (b) -6, (c) -5, (d) -4, (e) -3, (f) -2, (g) -1, and (h) 0. The meaning of symbols are 
as follows: “compr”; compressional heating, “H 2 cont”; cooling by H 2 collision-induced emission, 
“grain”; cooling by dust thermal emission, “H 2 form (diss)”; H 2 formation heating (dissociation 
cooling, respectively), which is indicated by “f”(“d”, respectively) in the Figure. Others are line- 
emission-cooling rates by indicated coolants. In particular, “[C&O]” means cooling by fine-structure 
lines of C’*', C, and O. In panel (d), “OH&H 2 O” is cooling by OH for nn < lO^crn”^ and by H 2 O 
for higher densities. Note that the compression and the H 2 formation are heating mechanism, while 
the other processes are cooling mechanism. 
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Fig. 3.— Evolution of H 2 fraction /hj = 2 y(H2) for the cases shown in Figure 1. Solid (dashed) 
lines indicate those for [Z/H]=-oo (Z=0), -5, -3, and -1 ([Z/H]=-6, -4, -2, and 0, respectively). Note 
that /h 2 = 1 means that hydrogen is fully molecular. 
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Fig. 4.— Evolution of HD fraction /hd = y(HD)/yD in low-metallicity clouds. Solid (dashed) lines 
show the cases of [Z/H]=- oo (Z=0), -5, and -3 ([Z/H]=-6, and -4, respectively). Note that /d = 1 
means that deuterium is fully molecular. 
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Fig. 5.— Important chemical network of (a) oxygen (b) carbon. All the C and O related reactions 
in the reduced network (see text) are shown. The indicated numbers are the same as Z reactions 
in Table 1, i.e., i is for reaction Zi. 
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Fig. 6.— Transformation of oxygen in metal-deficient clonds with metallicity (a)[Z/H]=-5, (b) -4, 
and (c) -3 in the fiducial cases. Concentration of species O, OH, H 2 O, and O 2 are shown against 
H number density nn. 
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Fig. 7.— Same as Figure 1, but with external radiation of CMB at l+z=20, i.e., radiation tem¬ 
perature Trad = 54.6K. Those of the fiducial cases are shown by dotted lines for comparison. 
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Fig. 8.— Same as Figure 1, but without gas-phase metals. The amount of dust is equal to the 
fiducial case with the same value of [Z/H] indicated. Also shown by dotted lines are those of the 
fiducial cases. Note that the total amount of metals are smaller than in the fiducial case with the 
same [Z/H] value because of the lack of gas-phase metals. 
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Fig. 9.— Same as Figure 1, but without dust grains. For clarity, fewer cases are shown here than 
in Figure 1. An amount of metals in gas phase is the same as the fiducial cases. Also shown by 
dotted lines are those of the fiducial cases. Note that the total amount of metals are smaller than 
in the fiducial case with the same [Z/H] value because of the lack of metals in dust. 
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Fig. 10.— Temperature evolution calculated by the minimum model, which includes fine-structure- 
line cooling and dust processes only, in addition to the metal-free gas microphysics. Among metal 
reactions, only the recombination of C (reaction Z28 in Table 1) are included. Also shown by dotted 
lines are those by the full network (the fiducial cases). 
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Fig. 11.— Temperature evolution with the reduced chemical network. The chemical reactions 
shown in Table 1 are included. Although those with full network are shown by dotted lines, all the 
lines are overlapping and cannot be recognized. 
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Fig. 12.— Abundances of important species by the reduced chemical network. Those with full 
network are shown by dotted lines. 
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Fig. 13.— Temperature evolution calculated by the model that includes only (a)[CII] and (b)[OI] 
cooling in addition to the metal-free processes. All C or O are assumed to be in the gas phase 
as in Bromm & Loeb (2003). The abundance of CII and 01 are assumed to be constant in time: 
the metal chemistry is not solved. Also shown by dotted lines are the fiducial cases with [Z/H] = 
—oo, 5,4, 3,2. 











Fig. 14.— Power of growth of the non-spherical deformation mode {m = 2). The amplitude of 
deformation £ evolves as oc p'^. The real part of exponent n is plotted as a function of effective 
adiabatic index 7 . 




1 



Fig. 15.— Evolution of elongation £ during the run-away collapse phase in the linear theory. If the 
elongation exceeds a critical value, the cores are assumed to fragment and the elongation is reduced 
as in the text. This effect is not included in this Figure. We assume that the run-away collapse 
begins after the first fragmentation, i.e., when the temperature begins to increase after the initial 
line-cooling phase. 
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Fig. 16.— The initial mass function of fragments for different metallicities. The number fraction in 
each bin is shown by the histogram. Each bin has width Alognn = 1 and centered where values of 
logrufrag are integer. The initial distribution of Ti/Tnl is (a)gaussian centered at 0.5 with variance 
IT = 0.1 and (b) uniform between 0 and 1. 





